Robust Low-Rank LCMV Beamforming Algorithms 
Based on Joint Iterative Optimization Strategies 

Rodrigo C. de Lamare 



o 



> 

m 
m 

O 

m 



X 



Abstract — This chapter presents reduced-rank hnearly con- 
strained minimum variance (LCMV) algorithms based on the 
concept of joint iterative optimization of parameters. The pro- 
posed reduced-rank scheme is based on a constrained robust 
joint iterative optimization (RJIO) of parameters according 
to the minimum variance criterion. The robust optimization 
procedure adjusts the parameters of a rank-reduction matrix, 
a reduced-rank beamformer and the diagonal loading in an 
alternating manner. LCMV expressions are developed for the 
design of the rank-reduction matrix and the reduced-rank beam- 
former. Stochastic gradient and recursive least-squares adaptive 
algorithms are then devised for an efficient implementation 
of the RJIO robust beamforming technique. Simulations for a 
application in the presence of uncertainties show that the RJIO 
scheme and algorithms outperform in convergence and tracking 
performances existing algorithms while requiring a comparable 
complexity. 

Index Terms — Adaptive beamforming, constrained optimiza- 
tion, robust techniques, reduced-rank methods, iterative methods. 



I. Introduction 

In the last decade, techniques have attracted a significant 
interest from researchers and engineers, and found applica- 
tions in radar, sonar, wireless communications and seismology 
m, The optimal linearly constrained minimum variance 
(LCMV) beamformer is designed in such a way that it min- 
imizes the array output power while maintaining a constant 
response in the direction of a signal of interest (Sol) O], 
El, El- However, this technique requires the computation 
of the inverse of the input data covariance matrix and the 
knowledge of the array steering vector Adaptive versions 
of the LCMV beamformer were subsequently reported with 
stochastic gradient (SG) ID, fSl and recursive least squares 
(RLS) Q algorithms. A key problem with techniques is the 
impact of uncertainties which can result in a considerable per- 
formance degradation. These mismatches are caused by local 
scattering, imperfectly calibrated arrays, insufficient training 
and imprecisely known wave field propagation conditions ^ . 

In the last decades a number of robust approaches have been 
reported that address this problem [8 |-[31 1. These techniques 
can be classified according to the approach adopted to deal 
with the mismatches: techniques based on diagonal loading 
ID, GOl, [131, 1 14 1, methods that estimate the mismatch or 
equivalently the actual steering vector ifTTI . lfT2l . ifTSl . and 
techniques that exploit properties such as the constant modulus 
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of the signals fT6l, flTl, fW\ and the low-rank of the in- 
terference subspace 19|,|19|-131 1. Furthermore, beamforming 
algorithms usually have a trade-off between performance and 
computational complexity which depends on the designer's 
choice of the adaptation algorithm [3 1, |6|. A number of robust 
designs can be cast as optimization problems which end up 
in the so-called second-order cone (SOC) program, which can 
be solved with interior point methods and have a computa- 
tional cost that is super cubic in the number of parameters 
of the beamformer This poses a problem for beamforming 
systems that have a large number of parameter and operate in 
time-varying scenarios, which require the beamformer to be 
recomputed periodically. 

A robust technique for short-data record scenarios is 
reduced-rank signal processing |fT9l - ll3n . which is very well 
suited for systems with a large number of parameters. These 
algorithms are robust against short data records, have the 
ability to exploit the low-rank nature of the signals encountered 
in beamforming applications and can resist moderate steering 
vector mismatches. These methods include the computation- 
ally expensive eigen-decomposition techniques |fT9l - ll20l to 
alternative approaches such as the auxiliary-vector filter (AVF) 
[21 J, [26,1, the multistage Wiener filter (MSWF) ll22l. Il24l. Il2l 
which are based on the Krylov subspace, and joint iterative 
optimization (JIO) approaches |23, lES), ED, EQI, ED- 
The JIO techniques reported in ll27l . Il28l . Il30l outperform 
the eigen-decomposition- and Krylov-based methods and are 
amenable to efficient adaptive implementations. However, ro- 
bust versions of JIO methods have not been considered so far 

In this chapter, robust LCMV reduced-rank beamforming 
algorithms based on constrained robust joint iterative opti- 
mization (RJIO) of parameters are developed. The basic idea 
of the RJIO approach is to design a bank of robust adaptive 
beamformers which is responsible for performing dimension- 
ality reduction, whereas the robust reduced-rank beamformer 
effectively forms the beam in the direction of the Sol and 
takes into account the uncertainty. Robust LCMV expressions 
for the design of the rank reduction matrix and the reduced- 
rank beamformer are proposed that can appropriately deal with 
array steering vector mismatches. SG and RLS algorithms 
for efficiently implementing the method are then devised. 
An automatic rank adaptation algorithm for determining the 
most adequate rank for the RJIO algorithms is described. A 
simulation study of the proposed RJIO algorithms and existing 
techniques is considered. 

This chapter is organized as follows. The system and signals 
models are described in Section II. The full-rank and the 
reduced-rank LCMV filtering problems are formulated in 
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Section III. Section IV is dedicated to the RJIO method, 
whereas Section V is devoted to the derivation of the adaptive 
SG and RLS algorithms, the analysis of the computational 
complexity, and the rank adaptation technique. Section VI 
presents and discusses the simulation results and Section VII 
gives the concluding remarks. 

II. System Model 

Let us consider a sensor-array system equipped with a 
uniform linear array (ULA) of M elements, as shown in Fig. 
1. Assuming that the sources are in the far field of the array, 
the signals of K narrowband sources impinge on the array 
{K < M) with unknown directions of arrival (DOA) 9i for 
/ = l,2,...,if. 




Output x{i) 



Fig. 1 . Block diagram of a sensor-array array system with interfering signals 

The input data from the antenna array can be organized in 
an M X 1 vector expressed by 



where 



r{i) = A{e)s{i) + n{i) 

A{e) = [a{e^),...,a{eK)] 



(1) 



is the M X K matrix of signal steering vectors. The M x 1 
signal steering vector is defined as 

T 



— 27rj cos (. 



, ,e 



-2Trj(M-l)j^ cost 



(2) 



for a signal impinging at angle 6i, I ^ 1,2, K, where dg — 
Ac/2 is the inter-element spacing, Ac is the wavelength and 
(.)^ denotes the transpose operation. The vector n{i) denotes 
the complex vector of sensor noise, which is assumed to be 
zero-mean and Gaussian with covariance matrix a^I. 

III. Problem Statement and Design of Adaptive 
Beamformers 

In this section, the problem of designing algorithms against 
steering vector mismatches is stated. The design of robust 
full-rank and reduced-rank LCMV beamformers is introduced 



along with the modelling of steering vector mismatches. The 
presumed array steering vector for the k-th desired signal is 
given by ap{9k) = a{9k) + e, where e is the M xl mismatch 
vector and a{6k) is the actual array steering vector which is 
unknown by the system. By using the presumed array steering 
vector ap{6k), the performance of a conventional LCMV 
beamformer can be degraded significantly. The problem of 
interest is how to design a beamformer that can deal with 
the mismatch and minimize the performance loss due to the 
uncertainty. 

A. Adaptive LCMV Beamformers 

In order to perform beamforming with a full-rank LCMV 
beamformer, we linearly combine the data vector r{i) — 



' M 



w 



] with the full-rank beamformer w 
to yield 



x{i) = w"r{i) (3) 

The optimal LCMV beamformer is described by the M x 1 
vector w, which is designed to solve the following optimiza- 
tion problem 

minimize i?[|it)'^r(i)p] = Rw 
subject to w^a{6k) = 1 
The solution to the problem in dU is given by O, ||4l 



(4) 



w. 



opt 



(5) 



afi{ek)R-'a{9k)) ' 

where a{9k) is the steering vector of the Sol, r{i) is the 
received data, the covariance matrix of r{i) is described by 
R = E[r{i)r^ (i)], {■)^ denotes Hermitian transpose and 
E[-] stands for expected value. The beamformer w{i) can be 
estimated via SG or RLS algorithms [3J. However, the laws 
that govern their convergence and tracking behaviors imply 
that they depend on M and on the eigenvalue spread of R. 

A reduced-rank algorithm must extract the most important 
features of the processed data by performing dimensionality 
reduction. This mapping is carried out by a M x D rank- 
reduction matrix Sd on the received data as given by 

f(z) = 5gr(i) (6) 

where, in what follows, all D-dimensional quantities are 
denoted with a "bar". The resulting projected received vector 
r{i) is the input to a beamformer represented by the D vector 
w = [iDi W2 ■ ■ ■ WjjY . The filter output is 

x{i) = w"r{i) (7) 

In order to design a reduced-rank beamformer w we consider 
the following optimization problem 

minimize i? [110^7^(1)1^] — ib^ Rib 

subject to w^d{dk) = 1 

The solution to the above problem is 



(8) 



R 



'Jopt 



-.Hi 



%)R 'diOk) 



{S%RSD)-^S%a{ek) 
aHSD{ek){SlRSD)-^Sla{ek) ' 



(9) 
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where the reduced-rank covariance matrix is ^ = 
E[f{i)f^{i)] = S^RSd and the reduced-rank steering 
vector is a{6k) = S^a{6k). The above development shows 
that the choice of Sd to perform dimensionahty reduction 
on r{i) is very important, and can lead to an improved con- 
vergence and tracking performance over the full-rank beam- 
former A key problem with the full-rank and the reduced- 
rank beamformers described in (|5} and (|9), respectively, is 
that their performance is deteriorated when they employ the 
presumed array steering vector ap{9k)- In these situations it is 
fundamental to employ a robust technique that can mitigate the 
effects of the mismatches between the actual and the presumed 
steering vector 



B. Robust Adaptive LCMV Beamformers 

An effective technique for is the use of diagonal loading 
strategies (S], lHO], [US, fU]. In what follows, robust full-rank 
and reduced-rank LCMV beamforming designs are described. 
A general approach based on diagonal loading is employed 
for both full-rank and reduced-rank designs. 

A robust full-rank LCMV beamformer represented by an 
M X 1 vector w can be designed by solving the following 
optimization problem 



minimize E[\w^r{i)\'^ 

subject to a{dk 



e'^WwW'^ =w"Rw + e^w"w 



(10) 



where is constant that needs to be chosen by the designer. 
The solution to the problem in (fTOl l is given by 



w. 



opt 



a"iek)iR + e^lM)-^api9k)) 



(11) 



where ap{dk) is the presumed steering vector of the Sol and 
Id is an A/-dimensional identity matrix. It turns out that 
the adjustment of needs to be obtained numerically by an 
optimization algorithm. 

In order to design a robust reduced-rank LCMV beamformer 
w we follow a similar approach to the full-rank case and 
consider the following optimization problem 



minimize E[\w" S^r{i)\^] +e^\\SDwf =w"S^RSdw 



IV. Robust Reduced-Rank Beamforming Based on 
Joint Iterative Optimization of Parameters 

In this section, the principles of the robust reduced-rank 
beamforming scheme based on joint iterative optimization of 
parameters, termed RJIO, is introduced. The RJIO scheme, 
depicted in Fig. 2, employs a rank-reduction matrix Soii) 
with dimensions M x D to perform dimensionality reduction 
on a data vector r{i) with dimensions M x 1. The reduced- 
rank beamformer w{i) with dimensions D x 1 processes the 
reduced-rank data vector r(i) in order to yield a scalar estimate 
x{i). The rank-reduction matrix S£)(i) and the reduced-rank 
beamformer w{i) are jointly optimized in the RJIO scheme 
according to the MV criterion subject to a robust constraint 
that ensures that the beamforming algorithm is robust against 
steering vector mismatches and short data records. 



r(i) 
M X 1 





r(z) 




Dimensionality 
Reduction 


Reduced-Rank 
Beamforming 
w(i) 




D X 1 


^ ' M X D 




- £> X 1 



x(i) 



Robust 
Algorithm 



Fig. 2. Block diagram of the RJIO scheme 

In order to describe the RJIO method, let us first consider 
the structure of the M x D rank-reduction matrix 



Snii) = [ si{i) I S2[i) \ ... \sD{i) 



(14) 



where the columns Sd{i) for d—\, . . . , D constitute a bank 
of D robust beamformers with dimensions M x 1 as given by 

Sdii) = [si,<i(i) S2,d{i) . . . smA'^)]'^ 

The output x{i) of the RJIO scheme can be expressed as a 
function of the input vector r(i), the matrix Soii) and the 
reduced-rank beamformer w(i): 



x{i) = w"{i)S^{i)r{i) = w"{i)f{i) 



(15) 



subject to S^ap{9k) 



1, 



The solution to the above problem is 



J opt 



(12) 



(13) 



where the tuning of requires an algorithmic approach as 
there is no closed-form solution and In is a £)-dimensional 
identity matrix. 



It is interesting to note that for D = 1, the RJIO scheme 
becomes a robust full-rank LCMV beamforming scheme with 
an addition weight parameter wb that provides an amplitude 
gain. For D > 1, the signal processing tasks are changed and 
the robust full-rank LCMV beamformers compute a subspace 
projection and the reduced-rank beamformer provides a unity 
gain in the direction of the Sol. This rationale is fundamental 
to the exploitation of the low-rank nature of signals in typical 
beamforming scenarios. 

The robust LCMV expressions for Soii) and w{i) can be 
computed via the following optimization problem 

minimize E[\w" {i)S^{i)r{i)\'^] + e^\\S D{i)w{i)\\^ = 

w"{i)S^{i)RSD{i)w{i) + e^w"{i)S^{i)SD{i)w{i) 

subject to w"{i)S^{i)ap{9k) = 1 

(16) 
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In order to solve the above problem, we resort to the method 
of Lagrange multipliers (3] and transform the into an uncon- 
strained one expressed by the Lagrangian 

^SD{^l w{^),e\^)) = E[\w^ {t)S^{t)r{t)\'] 

+ e^i^)w"i^)S'^i^)SD{^)wi^)) 
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(17) 



where A is a scalar Lagrange multiplier, * denotes complex 
conjugate. By fixing w{i), minimizing (fTTj l with respect to 
S]j{i) and solving for A, we get 



SD{^) 



{R + e\i)lM)-^ap{ek)w''ii)R^ 
w"{i)lQ^w{i)a^[ek){R{i) + e-^{i)lM)-^a.p 



(18) 



where R E[r{i)r" {%)] and Ra = E[w(i)w" (i)]. By 
fixing Soii), minimizing ([TtT i with respect to w{i) and 
solving for A, we arrive at the expression 



w{i) 



(Rji) + e^{i)S^{i)lDSDii))-^dp{ef. 
aH{ek){R{i) + e^{i)Sl{i)lDSD{i))-^a 



(19) 

where R{i) = £;[5'g(i)r(i)r^(j)5'D(i)] = E[r{i)f^{i)], 
<ip{0k) = S^{i)ap{9k)- Note that the filter expressions in 
(fTsT i and (fT9T l are not closed-form solutions for w{i) and 
Soii) since dTSI l is a function of w{i) and (fT9] l depends 
on Soii)- Thus, it is necessary to iterate (fTsl l and ([19) with 
initial values to obtain a solution. The key strategy lies in 
the robust joint optimization of the beamformers. The rank 
D and the diagonal loading parameter e^{i) must be adjusted 
by the designer to ensure appropriate performance or can be 
estimated via another algorithm. In the next section, iterative 
solutions via adaptive algorithms are sought for the robust 
computation of Soii), ^{i), the diagonal loading e{i) and 
the rank adaptation. 

V. Adaptive Algorithms 

In this section, adaptive SG and RLS versions of the RJIO 
scheme are developed for an efficient implementation. The 
important issue of determining the rank of the scheme with 
an adaptation technique is considered. The computational com- 
plexity in arithmetic operations of the RJIO-based algorithms 
is then detailed. 



A. Stochastic Gradient Algorithm 

In this part, we present a low-complexity SG adaptive 
reduced-rank algorithm for an efficient implementation of the 
RJIO method. The basic idea is to employ an alternating 
optimization strategy to update w{i) and the diagonal 

loading e^(i) By computing the instantaneous gradient terms 
of ([TtI i with respect to Soii), 'w{i) and e^(i), we obtain 

x*{i)r{i)w"{i) + e^{i)SD{i)w{i)w"{i) 



By introducing the positive step sizes /is, and /i^, using 
the gradient rules So^i + 1) = Soii) ~ [^s^ I^mv s*^{i)^ 
w{i + 1) = w{i) - fJ-w^jCMViv-ii) and e{i + 1) = e{i) - 
/i^ V£AfVe(i), enforcing the constraint and solving the result- 
ing equations, we obtain 

SD{^ + 1) = Sd{^) - A*. [x*{i)r{i)w"{i) + e{i)S Dit)w{i)w" (t) 

^ {a^{dk)ap{dk)y'ap{dk)w"{i){x*{t)a^{ek)r{i) + e{t) 

(21) 

wii + 1) = w{i) - ^l^{x*{i)S^i^)r{^) + e{i)S^ii)SD{i)w{t) 

+ {a^{9k)ap{0k))-'{x*{^)r"{^)SD{^)S'^{^)ap{0k) 

(22) 



- e{i)u 



e{i + 1) = e{i) - ii,w"{i)S'^(i)SD{i)w{i) 



(23) 



where x{i) = {i)S^{i)r{i). The RJIO scheme trades- 
off a full-rank beamformer against one rank-reduction matrix 
Snii), one reduced-rank beamformer w{i) and one adaptive 
loading recursion operating in an alternating fashion and 
exchanging information. 



B. Recursive Least Squares Algorithms 

Here, an RLS algorithm is devised for an efficient imple- 
mentation of the RJIO method. To this end, let us first consider 
the Lagrangian 



1=1 



(24) 



where a is the forgetting factor chosen as a positive constant 
close to, but less than 1. 

Fixing w{i), computing the gradient of (l24l i with respect to 
Soii), equating the gradient terms to zero and solving for A, 
we obtain 



P{t)ap{0k)a^{9k)SD{i-l) 
a^{0k)P{i)ap{0k) 



(25) 



where we defined the inverse covariance matrix P{i) — 
R~^{i) for convenience of presentation. Employing the matrix 
inversion lemma li3J, we obtain 



^P{i - l)r{i) 



W"{i)P{i - l)r(i) 



(26) 



P{i) = a-^P{i^l)-a-'^k{i)r"{i)P{i-l)+e'^{i)lM (27) 

where k{i) is the M x 1 Kalman gain vector We set P{0) = 
51 M to start the recursion of (l27t . where 5 is a positive 
constant. 

Assuming Soii) is known and taking the gradient of i 

solving for A, we obtain the D xl reduced-rank beamformer 

2\*S%{i)Kp{0k) 

P{i)SU^ap{ 



(20) 



w{i) 



a"{0k)SD{i)P{i)Sl{i)ap{0k) 



(28) 
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where P{i) = R ^(i) and R{i) = ELi "'"'^(O^'^CO is 
the reduced-rank input covariance matrix. In order to estimate 
P{i), we use the matrix inversion lemma [3J as follows 

a-^P{i - l)f{i) 



(29) 



I + a-^f" {i)P{i - \)f{i) 
P{i) ^ a~^P{i~l)-a-^k{i)f"{i)P{i^l)+e^{i)lD (30) 

where k{i) is the D xl reduced-rank gain vector and P{i) — 
R (i) is referred to as the reduced-rank inverse covariance 
matrix. Hence, the covariance matrix inversion R (i) is 
replaced at each step by the recursive processes (|29] and (l30T l 
for reducing the complexity. The recursion of (|30t is initialized 
by choosing P{0) = 61 jj, where 6 is a positive constant. The 
last recursion adjusts the diagonal loading according to the 
following update equation 



eii + 1) - e(^) - ^i.w'' {i)S^{i)S D{^)w{^) 



(31) 



The RJIO-RLS algorithm trades-off a full-rank beamformer 
with M coefficients against one matrix recursion to compute 
S£){i), given in (|25]l-(|27]i. one D x 1 reduced-rank adaptive 
beamformer w{i), given in (|28]|-(|30]|, and one recursion to 
adjust the diagonal loading described in (IJTT i in an alternating 
manner and exchanging information. 

C. Complexity of RJIO Algorithms 

Here, we evaluate the computational complexity of the RJIO 
and analyzed LCMV algorithms. The complexity expressed in 
terms of additions and multiplications is depicted in Table I. 
We can verify that the RJIO-SG algorithm has a complexity 
that grows linearly with DAI, which is about D times higher 
than the full-rank LCMV-SG algorithm and significantly lower 
than the remaining techniques. If D << M (as we will 
see later) then the additional complexity can be acceptable 
provided the gains in performance justify them. In the case 
of the RJIO-RLS algorithm the complexity is quadratic with 
AP and D^. This corresponds to a complexity slightly higher 
than the one observed for the full-rank LCMV-RLS algorithm, 
provided D is significantly less than M, and lower than the 
algorithms WC-SOC LIOJ and WC-ME UJJ. 

TABLE I 

Computational complexity of LCMV algorithms 



Algorithm 
LCMV-SG [41" 



LCMV-RLS tli 
RJIO-SG 

RJIO-RLS 

SMI PI 



Additions 



Multiplications 

:W + '2 



3M^ - 2M + 3 6M^ +2M + 2 



3DM + 4M 

+2D - 2 



SM^ - 
+3D^ 



M + 3 



5DM + 2M 
+5D + 2 

+ 3M 



2/3M^ + 3M2 2/3M^ + 5M^ 



In order to illustrate the main trends in what concerns 
the complexity of the proposed and analyzed algorithms, we 
show in Fig. 3 the complexity in terms of additions and 



multiplications versus the number of input samples M. The 
curves indicate that the RJIO-RLS algorithm has a complexity 
lower than the WC-ME U_l| and the WC-SOC LIOI, whereas 
it remains at the same level of the fuU-rank LCMV-RLS 
algorithm. The RJIO-SG algorithm has a complexity that is 
situated between the full-rank LCMV-RLS and the full-rank 
LCMV-SG algorithms. 



a) D=6 



b) D=6 




< — Full-rank-RLS 
RJIO-SG 
RJIO-RLS 
SMI 

e — WC-ME [ir 
.^WC-SOC [10] 



20 40 60 80 100 120 
M 



20 40 60 80 100 120 
M 



Fig. 3. 

M 



Computational complexity in terms of arithmetic operations against 



D. Rank Adaptation 

The performance of the algorithms described in the previous 
subsections depends on the rank D. This motivates the devel- 
opment of methods to automatically adjust D on the basis of 
the cost function. Differently from existing methods for rank 
adaptation which use MSWF-based algorithms |24| or AVF- 
based recursions [261, we focus on an approach that jointly 
determines D based on the LS criterion computed by the filters 
Soil) and WD{i), where the subscript D denotes the rank 
used for the adaptation. In particular, we present a method for 
automatically selecting the ranks of the algorithms based on 
the exponentially weighted a posteriori least-squares type cost 
function described by 

i 

C{Sd{i-1),wd{i-1)) = ^a*->g(^-l)5z,(^-l)r(^)|^ 

1=1 

(32) 

where a is the forgetting factor and W£)(i — 1) is the beam- 
former with rank D. For each time interval i, we can select the 
rank Dopt which minimizes C{SD{i ~ 1), ■u'd(* — 1)) and the 
exponential weighting factor a is required as the optimal rank 
varies as a function of the data record. The key quantities to be 
updated are the rank-reduction matrix Soii), the beamformer 
w !){{), the associated presumed reduced-rank steering vector 
CLpiOk) and the inverse of the reduced-rank covariance matrix 
P{i) (for the RJIO-RLS algorithm). To this end, we define 
the following extended rank-reduction matrix Soii) and the 
extended reduced-rank beamformer weight vector woii) as 
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follows: 



Snii) 



where Rs is the covariance matrix of the desired signal and 
Rj is the (Tovag^nce' 



Sl,l Sl,2 



SM.1 Sm,2 



Si, A, 



sm.d„ 



ani 



(33) 

The extended rank-reduction matrix Suii) and the extended 
reduced-rank beamformer weight vector wd{i) are updated 
along with the associated quantities a{6k) and P{i) (only for 
the RLS) for the maximum allowed rank Dmax and then the 
rank adaptation algorithm determines the rank that is best for 
each time instant i using the cost function in ( |32] |. The rank 
adaptation algorithm is then given by 



-Dopt = arg min 

Dmi„<d<r>„ 



C{SD{^~l),WD{^-l)) (34) 



where d is an integer, Z^min and Dmax are the minimum and 
maximum ranks allowed for the reduced-rank beamformer, 
respectively. Note that a smaller rank may provide faster 
adaptation during the initial stages of the estimation proce- 
dure and a greater rank usually yields a better steady-state 
performance. Our studies reveal that the range for which the 
rank D of the proposed algorithms have a positive impact 
on the performance of the algorithms is limited, being from 



3 to D„ 



for the reduced-rank beamformer 



recursions. These values are rather insensitive to the system 
load (number of users), to the number of array elements and 
work very well for all scenarios and algorithms examined. 
The additional complexity of the proposed rank adaptation 
algorithm is that it requires the update of all involved quantities 
with the maximum allowed rank Dmax and the computation 
of the cost function in (|32] i. This procedure can significantly 
improve the convergence performance and can be relaxed (the 
rank can be made fixed) once the algorithm reaches steady 
state. Choosing an inadequate rank for adaptation may lead 
to performance degradation, which gradually increases as the 
adaptation rank deviates from the optimal rank. 

VI. Simulations 

In this section, the performance of the RJIO and some 
existing beamforming algorithms is assessed using computer 
simulations. A sensor-array system with a ULA equipped 
with M sensor elements is considered for assessing the 
beamforming algorithms. In particular, the performance of 
the RJIO scheme with SG and RLS algorithms is compared 
with existing techniques, namely, the full-rank LCMV-SG |4] 
and LCMV-RLS [7], and the robust techniques reported in 
ifTOl . termed WC-SOC, and HIl, called Robust-ME, and the 
optimal linear beamformer that assumes the knowledge of 
the covariance matrix and the actual steering vector [2 |. In 
particular, the algorithms are compared in terms of the signal- 
to-interference-plus-noise ratio (SINR), which is defined for 
the schemes as 

w"{i)S'^(i)R,SD{i)w{i) 



SINR(i) 



w"{i)S%{i)RiSD(i)w(iy 



(35) 



matrix of the interference and noise 
in the envir jnmgj^t. Nate that for the full-rank schemes the 
i) = I M- For each scenario, 200 runs 
curves. In all simulations, the desired 
1, and the signal-to-noise ratio (SNR) 
-f. The beamformers are initialized 
0] and 5z,(0) = [II 

D) matrix with zeros in 



SINR(i) assumes. S"( 
-jirgjj^^^ Xo abtain the 
signal powe • iWS'^i^ 
is defined a 5 SNjl = 
as w{Q) 
where QoxM-d is a Z? x [M 



all experiments. 

In order to assess the performance of the RJIO and other ex- 
isting algorithms in the presence of uncertainties, we consider 
that the array steering vector is corrupted by local coherent 
scattering 



Or, 



= ) 



0: 



(36) 



fc=i 



where <^k is uniformly distributed between zero and 27r and 
Ok is uniformly distributed with a standard deviation of 2 
degrees with the assumed direction as the mean. The mismatch 
changes for every realization and is fixed over the snapshots of 
each simulation trial. In the first two experiments, we consider 
a scenario with 7 interferers at -60°, -45°, 30° -15°, 0°, 
45°, 60° with powers following a log-normal distribution with 
associated standard deviation 3 dB around the Sol's power 
level. The Sol impinges on the array at 30°. The parameters 
of the algorithms are optimized. 

N=32, K=8 users, SNR=15 dB 




- Full-Rank-SG 
— Full-Rank-RLS 
— WC-SOC [10] 
— WC-ME[11] 

- SMI 

- RJIO-SG 

- RJIO-RLS 

- Optimal 



6 8 10 

Number of snapshots 



12 



16 



Fig. 4. SINR performance of LCMV algorithms against rank (D) 
with M = 32, SNR = 15 dB, AT = 250 snapshots 



We first evaluate the SINR performance of the analyzed 
algorithms against the rank D using optimized parameters (/i^., 
and forgetting factors A) for all schemes and N = 250 
snapshots. The results in Fig. |4] indicate that the best rank 
for the RJIO scheme is _D — 4 (which will be used in the 
second scenario) and it is very close to the optimal full- 
rank LCMV beamformer that has knowledge about the actual 
steering vector An examination of systems with different sizes 
has shown that D is relatively invariant to the system size. 



7 



which brings considerable computational savings. In practice, 
the rank D can be adapted in order to obtain fast convergence 
and ensure good steady-state performance and tracking after 
convergence. 



N=32, SNR=15dB 




200 300 400 500 
Number of snapshots 



Fig. 5. SINR performance of robust LCMV algorithms against 
snapshots with M = 32, SNR = 15 dB 



We display another scenario in Fig. |5] where the robust 
adaptive LCMV beamformers are set to converge to the same 
level of SINR. The parameters used to obtain these curves are 
also shown. The curves show an excellent performance for 
the RJIO scheme which converges much faster than the full- 
rank-SG algorithm, and is also better than the more complex 
WC-SOC 1101 and Robust-ME [ill schemes. 



1=24, K=6 users, SNR=12dB 



M=24, K=6 users, SNR=12 dB 




- RJIO-SG-Auto 
RJI0-SG{D=3) 

- RJIO-SG{D: 



200 400 600 800 
Number of snapsfiots 



- RJIO-RLS-Auto 
RJI0-RLS{D=3) 
RJIO-RLS{D 



200 400 600 800 
Number of snapsfiots 



Fig. 6. SINR performance of RJIO-LCMV (a) SG and (b) RLS algorithms 
against snapshots with M = 24, SNR = 12 dB with rank adaptation 



In the next example, we consider the design of the RJIO-SG 
and RJIO-RLS algorithms equipped with the rank adaptation 
method described in Section V.D. We consider 5 interferers 
at -60°, -30", 0°, 45°, 60° with equal powers to the Sol, 
which impinges on the array at 15°. Specifically, we evaluate 



the rank adaptation algorithms against the use of fixed ranks, 
namely, D — 3 and D = 8 for both SG and RLS algorithms. 
The results show that the rank adaptation method is capable of 
ensuring an excellent trade-off between convergence speed and 
steady-state performance, as illustrated in Fig |6] In particular, 
the algorithm can achieve a significantly faster convergence 
performance than the scheme with fixed rank D — 8, whereas 
it attains the same steady state performance. 



N=24, K=6 users (1< i < 800), K=8 users (801 < i < 1600) SNR=12 dB 




Full-Rank-SG 

- Full-Rank-RLS 
-WC-SOC [10] 
-WC-ME[11] 

-SMI 

- RJIO-SG 

- RJIO-RLS 

- Optimal 



600 800 1000 
Number of snapsfiots 



1600 



Fig. 7. SINR performance of robust LCMV algorithms against the number 
of snapshots with M = 24, SNR = 12 dB in a non-stationary scenario 



In the last experiment, we consider a non-stationary scenario 
where the system has 6 users with equal power and the 
environment experiences a sudden change at time i — 800. 
The 5 interferers impinge on the ULA at —60°, —30°, 0°, 
45°, 60° with equal powers to the Sol, which impinges on the 
array at 15°. At time instant i = 800 we have 3 interferers 
with 5 dB above the Sol's power level entering the system with 
DoAs —45°, —15° and 30°, whereas one interferer with DoA 
45° and a power level equal to the Sol exits the system. The 
RJIO and other analyzed algorithms are equipped with rank 
adaptation techniques and have to adjust their parameters in 
order to suppress the interferers. We optimize the step sizes 
and the forgetting factors of aU the algorithms in order to 
ensure that they converge as fast as they can to the same 
value of SINR. The results of this experiment are depicted 
in Fig. |7] The curves show that the RJIO algorithms have a 
superior performance to the existing algorithms considered in 
this study. 

VII. Conclusions 

We have investigated robust reduced-rank LCMV beam- 
forming algorithms based on robust joint iterative optimization 
of beamformers. The RJIO reduced-rank scheme is based on a 
robust constrained joint iterative optimization of beamformers 
according to the minimum variance criterion. We derived 
robust LCMV expressions for the design of the rank-reduction 
matrix and the reduced-rank beamformer and developed SG 
and RLS adaptive algorithms for their efficient implementation 
along with a rank adaptation technique. The numerical results 
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for an application with a ULA have shown that the RJIO 
scheme and algorithms outperform in convergence, steady 
state and tracking the existing robust full-rank and reduced- 
rank algorithms at comparable complexity. The proposed 
algorithms can be extended to other array geometries and 
applications. 
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